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We address the presence of bound entanglement in strongly-interacting spin systems at thermal 
equilibrium. In particular, we consider thermal graph states composed of an arbitrary number of 
particles. We show that for a certain range of temperatures no entanglement can be extracted by 
means of local operations and classical communication, even though the system is still entangled. 
This is found by harnessing the independence of the entanglement in some bipartitions of such 
states with the system's size. Specific examples for one- and two-dimensional systems are given. 
Our results thus prove the existence of thermal bound entanglement in an arbitrary large spin system 
with finite-range local interactions. 
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I. INTRODUCTION 

The application of tools recently developed in the con- 
text of quantum information theory to problems histor- 
ically native from the area of many-body physics has 
helped us gain new insight about collective quantum phe- 
nomena [1 . In particular, the characterization of the 
entanglement properties of ground and thermal states 
of strongly-interacting spin Hamiltonians provides a de- 
scription of these systems from a novel and alternative 
perspective. Also, in general, it gives us information 
about how hard it is to simulate their dynamical and 
statical properties with classical resources [2 . In addi- 
tion, such characterization is crucial to determine when 
these systems can in turn assist as resources in some given 
quantum-information processing task. For all these rea- 
sons it is important to classify and characterize standard 
many-body models in terms of their entanglement prop- 
erties. 

However, the entanglement characterization in multi- 
particle systems turns out to be formidably hard. On 
one hand, the calculation of truly multipartite entangle- 
ment measures is in general extremely difficult, even for 
the pure-state case jUS]. On the other hand, for realis- 
tic many (and specially macroscopically many) -body sys- 
tems, the interaction with its surrounding environment 
can additionally never be neglected and mixed states 
have to be necessarily taken into account. Unfortunately 
though, the characterization of mixed-state entanglement 
is up to date very poorly developed even for the general 
bipartite case. One of the most frequent and important 
type of such interactions - the one on which we will fo- 
cus here - takes place when the system is embedded in 
a thermal bath at temperature T and reaches thermal 
equilibrium with it, a process called thermalization. This 
process typically causes the system to lose quantum co- 
herence, gain entropy and in most times also lose entan- 
glement. 

An important step forwards in the characterization of 
mixed-state entanglement was to recognize two different 
types of entanglement: distillable and bound entangle- 



ment [5 . An entangled state p is said to be distillable if 
it is possible, by means of local operations and classical 
communication (LOCC), to obtain from p (or, more pre- 
cisely, several copies thereof) pure-state entanglement. 
Entangled states for which this task is impossible are said 
to be hound- entangled. Bound-entangled states were for 
some time believed to be useless for quantum informa- 
tion processing. Nevertheless, they are nowadays known 
to be useful in some practical situations [6^. 

Historically, most examples of bound entangled states 
have been provided without relying on any operational 
recipe (see e.g. Ref. [4 ). Yet this kind of entanglement 
has been recently found to arise in natural processes, in 
particular thermalization [7j [H [TOl [11 and other dy- 
namical decoherence processes JJ. For thermal states, 
the presence of bound entanglement was found in chains 
of harmonic oscillators in the thermodynamic limit, that 
is for an arbitrary number of oscillators [7] . On the other 
hand, in the case of fermionic systems the existence of 
thermal bound entanglement has been shown only for 
small systems of up to 12 spins [7[|8l[9] or for models in- 
volving non-local interactions between an arbitrary num- 
ber of spins [10 . The main contribution of the present 
paper is thus to show the existence of bound entangled 
thermal states in strongly-correlated systems of arbitrary 
number N of spin-^ particles with local interactions, in 
particular in the limit N — > oo. 

We show the presence of thermal bound entanglement 
for the exemplary family of Hamiltonians through which 
graph states are defined [12]. These Hamiltonians are 
always frustration-free and typically, depending on the 
graph, also local (meaning that not all particles inter- 
act simultaneously) and of finite-range (meaning that no 
particle interacts with another one infinitely far away) 
[13]. Graph states constitute an extremely important 
family of states from practical and fundamental points 
of view. They include cluster states, which are resources 
for universal measurement-based quantum computation 
[14 , codeword states for quantum error correction [15], 
and the well-known GHZ states, that are resources for se- 
cure quantum communication [16]. Moreover, this fam- 
ily of states can be used in quantum non-locality tests 



2 



[T7l[l8l[ia[20]. 

The motivation of this work is thus twofold: From a 
practical point of view, to establish the range of temper- 
atures for which a thermal state of an interesting many- 
body model proves useful as a resource for some poten- 
tial quantum information processing task. And, from a 
fundamental viewpoint, to take this particular case as a 
concrete example within a broader investigation of the 
properties of quantum correlations of strongly-correlated 
systems undergoing open-system dynamics. 

II. THERMAL GRAPH STATES AND 
DEPHASED GRAPH STATES 

In this section we define the states under scrutiny and 
settle the notation. Let us then start by defining a math- 
ematical graph G = {V, £} as the union of the set V of 
vertices z G V with the set £ of edges {^, j} G £ connect- 
ing each vertex i to some other j, being 1 < j < N. 
Next, for each graph G we define Hamiltonian H acting 
on N spin- 1/2 particles as 

1 ^ 

H = --^B,X,^ Z,, (1) 

i=l jeAfi 

where Bi > are arbitrary (strictly positive) coupling 
strengths in arbitrary units, Xj^ and Zj. are the usual 
Pauli operators acting on particle /c, and A//c denotes 
all neighbouring particles of k - i.e. each particle whose 
graph representation in G is a vertex j directly connected 
to k by some edge {j^k} G £. Hamiltonian ([T]) involves 
m-body interactions, where m is given by the maximum 
connectivity of the graph G. Also, since each and all 
of the N terms in summation ([T]) commute with each 
other, the eigenstates - here denoted by \G^^ ... ^^), with 
/if = or 1, for 1 < i < N - of each local term are 
also eigenstates of the whole summation. This implies 
that Hamiltonian (IT]) is frustration- free, meaning that 
the ground state -\Go ... o) - minimizes the energy of 
each term in H. 

Ground state |Go ... o) is the unique (non-degenerate) 
ground state of Hamiltonian ([T]) with eigenenergy 
— ^"^iLiBi^ whereas states \G^^ ... hn) ~ related 
to the former by the local-unitary transformation 
... AXTv) = ®^=i^i^'\Go ... o) - are eigenstates of 
([T]) with Xli^i Mi excitations and associated eigenener- 
gies The eigenstates \G^, ... 

form a complete orthogonal basis of the A/'-qubit Hilbert 
space. Since they are local-unitarily related, they all 
possess exactly the same entanglement properties, ex- 
tensively studied in Ref. [T2j 121] and references therein. 
For historical reasons, the ground state |Go ... o) has been 
taken however as the defining state for the so-called graph 
state. Let us recall that there exists also an alternative 
operational definition for such graph state: It can be 
physically produced initializing qubits in the super- 
position 1+) = (|0) + |1))/V^ and subsequently applying 



control-Z gates CZij = e*i(^^^^J"^^"^^+^) onto each 
pair of neighboring qubits defined by the graph G. Math- 
ematically: 

N N 

|Go... o)=(g)(g)CZ,,(g)|+)fc. (2) 

i=l jeAfi k=l 

We are now in condition to introduce the thermal 
graph state associated to G as the thermal state of Hamil- 
tonian ([T]): 

where T is the temperature of some bath with which 
our system of interest has reached thermal equilibrium 
(Boltzmann's constant is set as unit ks = I throughout). 

Alternatively, pr can be defined as a decohered graph 
state. In order to show that we need to introduce the 
completely-positive map A, acting on any A/'-qubit den- 
sity matrix p as 

A(p)=Pi0... P^(p), (4) 

as the composition of local, independent channels P^, 

Vi{p)={l-^)p+^Z,pZi, (5) 

with < Pi < 1. Local channel Vi describes the phys- 
ical process in which, with probability p^, an undesired 
TT-phase shift is experienced by qubit i and, with proba- 
bility I — Pi, the system is left untouched. Such process 
is present in situations where, with probability pi, there 
is complete loss of quantum coherence but without any 
population exchange. In the context of decoherence, map 
A in turn is often referred to as local (or individual) de- 
phasing (or phase damping). 

It was shown in Ref. [22 that thermal state ([3| - for 
the particular case of constant couplings Bi = B - can 
also be alternatively obtained by individually dephasing 
graph state (|2| with equal probabilities pi = p = ^_^^b/t • 
As shown in Appendix [A| this property also holds for 
general thermal states of Hamiltonian ([T]) with arbitrary 
couplings. Indeed, 

PT = A(|Go... o)(Go... ol), (6) 

with A defined according to Eqs. (|4| and ([5|, but with 
the additional constraint that the local dephasing prob- 
abilities satisfy 

The equivalence mathematically expressed in Eqs. (|6| 
and ([?]) establishes a very interesting connection between 
a collective decoherence process [thermalization of sys- 
tems governed by graph- state Hamiltonians as ([T])] and a 
local one [individual dephasing of systems initialized in 
graph states as ([2|]. 
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III. APPEARANCE OF MULTIPARTITE 
BOUND ENTANGLEMENT IN THERMAL 
GRAPH STATES 



As mentioned before, the characterization of entangle- 
ment in multiqubit systems is formidably hard even for 
pure states. However, a useful tool for the evaluation 
of the amount of entanglement contained in decohered 
graph states was developed in Ref. [23] for an important 
family of decoherence processes (see also Refs. [22 , 24 ). 
Given a certain multipartition of a graph G, the ma- 
chinery developed in [23 allows to map the calculation 
of the entanglement of a decohered (mixed) graph state 
to the average entanglement of several effective systems, 
constituted only by the so-called boundary qubits - the 
ones lying on the border of the multipartition and having 
neighbours on the other side of the border. Solving the 
former problem involves an optimization over a param- 
eter space exponentially large with whereas solving 
the latter involves only an optimization over the bound- 
ary qubits, a task that requires always exponentially less 
memory space and usually also considerably less compu- 
tational time, specially when N is large as in the ther- 
modynamical limit. 



For our case of interest, individual dephasing, this for- 
malism works even better as the entanglement contained 
in an arbritrary multipartition of a locally-dephased 
graph state is equivalent not to the average entangle- 
ment of an ensemble of smaller effective boundary sys- 
tems but just one. Furthermore, such effective system 
is simply composed by the locally-dephased original sys- 
tem itself but without all non-boundary qubits. The key 
point behind this idea is that all the control- Z gates that 
define |Go ... o) in Eq. ([2| commute with the dephas- 
ing map Q [25^. Hence, the order in which channels £i 
and the control- Z gates are applied on the product state 
(8)fc=i to obtain A(|G ... o){Go ... o|) is irrelevant. 
In particular, state ([6| is also obtained if the control-Z 
gates act after the dephasing channels. Now, all GZ's not 
crossing any boundary are local unitary operations with 
respect to the multipartition (see Figs. [l]and|3]for simple 
examples). Thus, because every entanglement quantifier 
is invariant under local unitary operations, as far as what 
concerns the amount of entanglement in the multiparti- 
tion one can simply forget about these non boundary- 
crossing GZ's. 



In what follows we apply this idea to some well-known, 
paradigmatic examples of graphs to show that there ex- 
ists a range of temperatures where the associated thermal 
graph state possesses multipartite bound entanglement. 
We do it first for the case of constant couplings Bi = B 
to transmit the essential idea clearly and then move to 
the arbitrary-coupling case. 



A. The linear cluster with equal couplings 

Let us start by the simplest example: the linear clus- 
ter state. Here the defyning graph G is the linear graph 
sketched in Fig. [l] We denote its thermal state as Pidt- 
First we consider a bipartition of the system into two 
contiguous blocks of spins (Fig. [l]A.), say from qubit 
i to the left (grey) and from qubit z + 1 to the right 
(white). We can easily see from the figure that all control- 
Z gates but one (in blue) act as local unitary opera- 
tions, and thus have no effect on the entanglement in 
the bipartition considered. As a result, the entire entan- 
glement between any two contiguous blocks of spins in 
Pidt is equivalent to that of the simple two-qubit ther- 
mal graph state in boundary pair i-i -\- 1. Then, by im- 
posing the entanglement between i and z + 1 to vanish 
we can establish the critical temperature for which the 
entire thermal cluster is separable with respect to the 
bipartition under scrutiny. Any entanglement quantifier 
valid for two-qubit mixed states would do for this aim, 
so we choose the simplest one to calculate: the negativ- 
ity [26]. The negativity Neg{p) of a state p is the sum 
of the absolute values of the negative eigenvalues of p^, 
where p^ is the partial transposition of p according to 
some bipartition. The negativity Neg^ij^i^piDr) ^^le 
thermal state of pair i\i ^ 1 is readily calculated to be 
Negi\ij^i{piDT) = i(2 - ^Pi - 2pi+i +PiPi+i). For the 
case of constant couplings Pi = p^+i = p = i_^^b/t condi- 
tion Negi\i_^i{piDT) = leads to the critical temperature 
i22j 



-B 



I.IB. 



(8) 



For T > T'^i\i-\-i{B)^ since the partition is separable, 
it is not possible to extract any type of entanglement be- 
tween any two contiguous blocks joined at spins i and i-\-l 
from a thermal ID graph state by applying contigous- 
block-local operations (arbitrary operations acting on 
spins 1 to z and on spins z + 1 to N) . This imposes strong 
restrictions on the operations one needs to apply to Pidt 
to distill some entanglement (if any) from some other of 
its multipartitions. For instance, if individual local oper- 
ations on each spin are applied, no entanglement can be 
distilled from this bipartition of Pidt T > T^i|i+i(5), 
for these are a particular case of contigous-block-local 
operations. Now, since here we have taken all coupling 
strengths equal, the critical temperature for separability 
of contiguous blocks is the same for all i. This means 
that for T > T'^^i_^i{B) all contiguous blocks of Pidt 
are in a separable state. So, no entanglement between 
any two particles can be extracted by LOCC, as for any 
two particles a contiguous-block bipartition can be found 
in which each particle lies on a different side of the par- 
tition and is therefore separable from the other. Ergo, 
no entanglement at all can be extracted from Pidt 
T>T-,|,+i(5) by LOCC. 

We now find another (non-contiguous) family of 
bipartitions for which the separability temperature. 



4 



A. 







i 






# — 


• 


• 


— • — • — 


• # 

J 


B. 


r 


i - 1 


r 

i 






#— 


• — 

) 




— • • — 

K 


• # 



FIG. 1: Two possible bipartitions of the linear cluster: the 
system is split into two subpartitions, represented by the grey 
and white regions. Note that all control-Z gates correspond- 
ing to black edges act locally with respect to the bipartition 
adopted, and thus do not change its entanglement. A. The 
systems is divided in two contiguous blocks of spins. The 
entanglement in this bipartition is equivalent to the entan- 
glement between the two qubits in blue (see text also). B. 
Another possible bipartition, this time non-contiguous blocks 
of spins are part of the same single subpartition (the white 
one). In this case all the entanglement is contained in the 
three-qubit boundary system shown in blue. 



i\i-\^i^\{B\ is strictly larger than T^i|i+i(5). This 
suffices to prove the non- separability oi P\bt~ there- 
fore the presence of multipartite bound entanglement 
in - for a range of temperatures T^i|i+i(5) < T < 



{B) 



Consider for instance the entangle- 



ment between the i-th qubit in the chain and all the 
rest (see Fig. [i]B). In this case, one can ignore all but 
two control-Z gates in the calculation of the entangle- 
ment. So, for these partitions, the entanglement (again 
quantified here by the negativity) of Pidt equiva- 
lent to that of the central particle vs. its two neigh- 
bours in a linear thermal cluster state of only three 
qubits. This negativity vanishes at a critical tempera- 
ture T^i\i_i^ij^i{B) that turns out to be strictly larger 
than T^,|,+i'(5) (T^,|,_i,,+i(5) > 1.65, see Fig. In 
the range T^i|i+i(5) < T < T^i|i_i,i+i(5), even tTiough 
Pidt possesses entanglement, we already know that none 
of it can be distilled through local operations assisted by 
classical communication. Such entanglement can only be 
extracted if particles i — 1 and i + 1 interact, which is of 
course not a local operation. Therefore, in this range of 
temperatures thermal cluster state Pidt possesses mul- 
tipartite bound entanglement. 

It is important to stress that none of the latter results 
or conclusions depends at all on the size N of the graph 
This allows us to guarantee that thermal bound entan- 
glement is also present in macroscopic specimens of these 
graphs. This independence of the entanglement on the 
graph's size is precisely the key point behind the method 
we used to simplify the calculation of the negativity of 
Pidt' general, in order to calculate the negativity of 




Thermal Bound 
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FIG. 2: Negativities of the thermal linear cluster state Pidt^ 
as a function of temperature T/B. The considered biparti- 
tions are those of Fig. 1: any two contiguous blocks, in red 
and any qubit versus the rest, in blue. None of these curves 
depends on the size N of the graph. Therefore such thermal 
bound entanglement is also present in the macroscopic ther- 
modynamical limit. Besides, we also display in dashed line 
the negativity of the even-odd partition, where particles with 
even label belong to subsystem A and particles with odd label 
to subsystem B. In this case, the negativity does depend on 
the system's system and the plot is done for N — 12. Our nu- 
merical investigations suggest that the even-odd partition is 
the most robust bipartition, i.e. the one with the highest crit- 
ical temperature of vanishing negativity. The shadded region 
shows the range of temperatures where the thermal system 
possesses bound entanglement (see text). 



an A^-qubit mixed state, one would need to diagonalize 
a 2^ X 2^ matrix, which already for a few tens of qubits 
cannot even be written down by a current classical com- 
puter. In the example studied in this subsection, this 
method has enabled us to obtain results for arbitrarily 
large systems calculating only negativities of two-qubit 
and three-qubit systems. 



B. The 2D square cluster with equal couplings 

Our next example is the 2D thermal cluster state P2Dt^ 
for which the associated graph G is a ^/N x \/N square 
lattice (see Fig. Isl). The columns of G are labeled by 
index 1 < i < and the rows by index 1 < j < a/]V. 
At T = this state is known to be a universal resource 
for one-way quantum computation. As a consequence, 
understanding the entanglement properties of this model 
under realistic noisy conditions is of course very impor- 
tant from a practical point of view. 

To prove the existence of thermal multipartite bound 
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FIG. 3: Two different bipartitions of a 2D cluster. A. The system is divided vertically into two parts. The only control-Z 
gates not acting as local unitary operations are the boundary-crossing ones, shown in blue. Thus the entanglement in any 
such bipartition of the 2D cluster in a thermal state is equivalent to that of several copies of a two-qubit thermal cluster 
state (represented in blue). B. The cluster is partitioned into a central spin and all other ones. The entanglement is this 
time equivalent to that between the central spin and only its four neighbours in a five-qubit thermal graph state in a star 
configuration (also in blue). 



entanglement in this example we follow exactly the same 
reasoning as in the previous subsection: establish a re- 
gion of temperatures where the system is separable with 
respect to any two contiguous blocks - and therefore 
nondistillable with respect to (contiguous-block) local op- 
erations while still being entangled in other partitions. In 
Fig. [3] A we consider a vertical partition of the system 
from column i to the left (grey) and from column i -\- 1 
to the right (white). Following the same steps as be- 
fore, the calculation of entanglement in this bipartition 
can be reduced to that of a product of the a/TV two- 
qubit thermal cluster states lying on the boundary. In 
this way we see that the critical temperature for sepa- 
rability (of all contiguous blocks) is again given by ([8|. 
Once again we consider another family of bipartitions, 
say, any qubit ij inside the lattice {1 < i < \/N and 
1 < j < ^/N) vs. all the rest of the qubits (see Fig. [3]b). 
The entanglement in this bipartition is equivalent to that 
between the qubit ij and its four neighbours, in a ther- 
mal graph state in a star configuration. The temperature 
for which the negativity vanishes in this case happens to 
be T^ij\i_i^ij^i^jj^i^j_i{B) ^ 2.55, which is again strictly 
larger than T^i|i+i(5). 



C. Unequal Hamiltonian couplings 



Consider now arbitrary couplings Bi in Hamiltonian 
([1]) with average value 5, that is: Yl!i=i = We 
first restrict ourselves to the familiar one-dimensional 
graph studied in subsection [III A[ as this already qual- 
itatively captures all the essential changes that can ap- 
pear when different couplings are present. The contigous- 
block separability condition Negi\i_^i{piDT) = ~ 



2pi - 2pi^i + PiPi^i) = now involves Pi = -^^^bjt 

= i^e^^+i/^ ' Thus the critical temperature is im- 
plicitly expressed by the equation 



1 



+ e T . (9) 

Equation ([9| above has a unique real solution T = 
^i+i)^ but is however non-invertible in gen- 
eral. Clearly, critical temperature T^i\ij^i{Bi^Bij^i) is 
symmetric under the exchange of Bi and Also, 
numerical inspections immediately show that it is a 
monotonously growing function of Bi (^i+i) for fixed 
Bij^i {Bi). Furthermore, under some contraints, as for 
example Bi + ^i+i constant, the critical temperature is 
maximal when the couplings are equal. 

As before, we compare T^i|i+i(5i, 5^+i) with the 
critical temperature T^^|^_i ,j+i(5i_i, 5^, 5^+1) corre- 
sponding to the bipartition which distinguishes qubit 
i versus the rest. For all cases we have studied, 

i\i_i^ij^i{Bi_i^ Bi^ ^i+i) turns out to be strictly larger 
than T''i\ij^i{Bi, Bi j^i) M I < i < N . If this is true in 
general, it implies that neither the appearance, nor the 
range of temperatures, of bound entanglement will be 
considerably affected by small deviations of Bi from the 
mean value B. 

However, the critical temperatures can be consider- 
ably sensitive to the values of the coupling constants at 
each site i when the deviations are large (say, of the or- 
der of B itself). In such case, it might as well happen 
that T^i\ij^i{Bi,Bij^i) < T^i\i_i^i^i{Bi_i, Bi, Bi^i) < 
T'^k\k+i{Bk,Bk+i) < T''/e|/c-i,/c+i(5fe-i,5/c,5fe+i), for 
some i ^ k. Remarkably, also in these cases the sys- 
tem displays a finite range of temperatures for which 
it is bound-entangled. To see this, it suffices to con- 
sider the site imax for which the critical temperature 
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of contigous-block separability is the largest of all, i 

y 1 < i < N. By the same reasonings as before, from 
on Pidt is non-distillable, but it is entangled up to 

T T j _ \ ^iqmax~\~^^ max 1 ' max ^ m a cc ~l~ 1 ) 



ax I ''max 



ix I 



max * 



Therefore (assuming again that 



is always true), 
-i,imax+^) state is 



for T e [T^ 

bound-entangled. The main conclusion to draw from the 
considerations in this paragraph is that rather than a pe- 
culiarity of the (ideal) case of equal coupling strengths, 
the presence of bound entanglement in a finite range of 
temperature appears to be a general phenomenon for (ar- 
bitrarily large) thermal graph states 

Finally, let us stress that similar results hold also for 
other graphs with different couplings. There again differ- 
ent partitions will give rise to different critical tempera- 
tures implying again the presence of bound entanglement 



along the lines of Sec. IIIB 



IV. CONCLUSION 



It is also worth mentioning that we used the negativity 
of some bipartitions of the system to detect inseparabil- 
ity of the thermal state. So our method can not detect 
bound entanglement in the case that all bipartitions of 
a system have positive partial transpositions. This was 
in turn found in the case of small spin systems [8 . We 
leave the existence of such kind of bound entanglement 
in macroscopic systems as an open problem. 

An interesting question concerns the utility of ther- 
mal graph states for information processing. For the re- 
gion of temperatures such that these states are distillable, 
one could in principle apply first a distillation protocol 
{e.g. see [22]) before using the state as a resource for 
quantum information. However, for the regions of bound 
entanglement, not even this experimentally-demanding 
strategy would work. 

Since in any practical implementation the temperature 
is always non-null, thermal bound entangled graph states 
are ideal probes to explore the limitations of realistic (ex- 
perimentally feasible) measurement-based quantum com- 
putation. This question is certainly of great interest and 
can be the subject of further analysis. 



Considering thermal graph states we have shown the 
presence of bound entanglement in systems containing a 
macroscopic number of spins with finite range interac- 
tions. This result extends previous results [7 for bosonic 
chains to fermionic systems. Our findings suggest that 
thermal bound entanglement could manifest also in more 
general systems, since this seems to be a robust feature 
against, in particular, to modifications in the symmetries 
of the Hamiltonian. 

In this paper we considered systems with three (or 
more) body interactions. Thus a natural question arises 
regarding the presence of macroscopic thermal bound en- 
tanglement also in spin models containing only 2-body 
interactions. This was indeed found in the case of har- 
monic oscillators systems [7 . Actually, we expect ther- 
mal bound entanglement to be a common phenomenon 
of general many-body systems since it is very unlikely 
that the negativities of all possible bipartitions of a sys- 
tem vanish at the same temperature. However, proving 
that for systems in the thermodynamical limit turns to 
be nontrivial. 
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APPENDIX A: THERMAL GRAPH STATE AS A 
DEPHASED GRAPH STATE FOR ARBITRARY 
COUPLINGS 

Thermal state Jsl) expressed in the eigenbasis 
of Hamiltonian Q, {\G^^ ... /i^)}, reads (dis- 
regarding its normalization) pT = e~^^^ = 

Z^/ii ... /iAr=0 ^ I^Mi ••• M7v/\'^Mi ••• MTV r 

Next we explicitly evaluate dephased state 
^{\Gq ... o)(Go ... o|) using Eqs. Q and ([s]) : 



A(|^o ... o)(Go ... ^\)=V^® ... D^(|Go ... o)(Go ... o|) 



1 N 
Ml ••• fJ'N=0 i=l 



where "| I2" stands for "modulo 2". The latter is 
equal to the former expression for pT if and only if 
each and all of the terms in the summation are equal. 
That is, if and only if n^^i (1 - Pi/2)l^^"^l^ (Pi/2)I^^|2 = 



g2T T^Z^iBii 1)'^^^ which when Eq. ^ holds -and using 
the fact that -^_^^b^/t = g-B^/2T^gBi/2T ^^nd 1 — -^_^^b^/t = 
g-B^/2T^gBi/2T lu turn reads H^^x g-Bi/2T^gB^/2T — 
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Q2T^7=i^i^ 1)'^* which can in turn be immediately tor. 
checked to be true up to a constant normahzation fac- 
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